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We show results for the T spectrum calculated in lattice QCD including for the first time vacuum 
polarization effects for light u and d quarks as well as s quarks. We use gluon field configurations 
generated by the MILC collaboration. The calculations compare the results for a variety of u 
and d quark masses, as well as making a comparison to quenched results (in which quark vacuum 
polarisation is ignored) and results with only u and d quarks. The b quarks in the T are treated 
in lattice Nonrelativistic QCD through NLO in an expansion in the velocity of the 6 quark. We 
concentrate on accurate results for orbital and radial splittings where we see clear agreement with 
experiment once u, d and s quark vacuum polarisation effects are included. This now allows a 
consistent determination of the parameters of QCD. We demonstrate this consistency through the 
agreement of the T and B spectrum using the same lattice bare b quark mass. A one-loop matching 
to continuum QCD gives a value for the b quark mass in full lattice QCD for the first time. We 
obtain m,^'^{m^'^) = 4.4(3) GeV. We are able to give physical results for the heavy quark potential 
parameters, ro — 0.469(7) fm and ri — 0.321(5) fm. Results for the fine structure in the spectrum 
and the T leptonic width are also presented. We predict the T — r/b splitting to be 61(14) MeV, 
the T' — rj'b splitting as 30(19) MeV and the splitting between the ht and the spin-average of the 
Xb states to be less than 6 MeV. Improvements to these calculations that will be made in the near 
future are discussed. 



I. INTRODUCTION 

The calculation of the spectrum of bottomonium and 
charmonium states is a key one for lattice QCD. for sev- 
eral reasons: 

• There are many 'gold-plated' (narrow, stable and 
experimentally well-characterised) states whose 
masses can be calculated accurately in lattice QCD. 

• Splittings in the spectrum have particularly good 
properties to use for determining the scale of the 
theory, Aqcd- In lattice QCD determining this 
scale is equivalent to knowledge of the lattice spac- 
ing and it is an important first calculation that 
must be done before quoting other dimensionful 
quantities on the lattice such as hadron masses, or 
giving a scale to as- 

• The bottomonium and charmonium calculations 
provide a test of the lattice QCD actions used for 
the heavy b and c quarks. These can then be used 
for the B and D physics results needed for the 
experimental programme charged with testing the 
consistency of CP violation in the Standard Model. 
Lattice QCD numbers are critical, for example, to 
push down the limits on determination of the ele- 
ments of the CKM matrix to the few percent level. 
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• The control of systematic errors from the lattice 
for heavyonium masses is well developed and the 
results are statistically very precise. Vacuum po- 
larisation effects for heavy quarks are not included 
because they are known perturbatively to have very 
little effect. This makes the tuning of the heavy 
quark mass much simpler than that of light quark 
masses. The heavyonium spectrum is sensitive to 
u, d, and s quark vacuum polarisation effects, but 
is not very sensitive to the actual values of the 
light quark masses, once they are low enough. All 
of these points make calculations of heavyonium 
masses in the presence of light quark vacuum polar- 
ization relatively simple, and a good test of lattice 
methods. 

We need to achieve statistical and systematic errors of 
a few percent in order for lattice QCD to provide use- 
ful information for the particle physics programme not 
obtainable by other methods. New results Q, 1^,1^1^,1^ 
indicate that this is starting to be possible, building on 
many years of work by the lattice community to under- 
stand the sources of error in lattice calculations. The 
developments that have enabled this progress are: im- 
proved actions, the inclusion of light sea quarks (i.e. light 
quark vacuum polarization effects) and effective theories 
for handling heavy quarks. In particular the inclusion of 
u, d, and s sea quarks with light u and d quark masses 
has produced agreement with experiment for a range of 
hadron masses for the first time A key component 
of this, for the reasons above, has been a set of results 
on the T spectrum and we describe these in more detail 
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here. 

Section |n] describes details of the lattice calculations 
and section imi gives the results obtained. In sections Hvl 
and ^ wc describe the analysis that allows the extrac- 
tion of the lattice spacing and wc compare the spectrum 
obtained to experiment. In section IVll we describe the 
determination of the b quark mass from the results. In 
sections I V III and I V IIII we discuss the fine structure in the 
spectrum and the determination of the leptonic width 
of the T and its radial excitations. We conclude in sec- 
tion IIXI with a discussion of further improvements that 
will be made, particularly to these last two calculations. 



lattice spacing values around 0.12fm and 0.09fm |3, and 
recently a further, coarser set, has become available with 
a lattice spacing around O.I7fm. 



The physical volumes of the configurations are large, 
(2.5fm)'^, so we do not expect significant errors from hav- 
ing a finite volume. Some analysis of finite volume errors 
has been done for light hadrons by the MILC collabo- 
ration 01 and no significant effect was found on these 
volumes. T states are smaller than light hadrons in gen- 
eral so we expect even less of an effect here. 



II. THE LATTICE CALCULATION 

The bottomonium spectrum is obtained by calculating 
correlators for bottomonium hadrons on configurations of 
gluon fields generated by a Monte Carlo procedure. The 
gluon field configurations in an ensemble arc generated 
with a probability distribution given by exp^—SgcD) 
where Sqcd is the QCD action. Averaging over results 
on an ensemble then gives the Feynman Path Integral 
result for the bottomonium correlator. First we describe 
the different ensembles of gluon field configurations that 
we have used and then the method for calculating the 
bottomonium correlators. We will then estimate the size 
of systematic errors in the calculation and finally discuss 
the fitting procedure which enables us to obtain energies 
and amplitudes for different states from the correlators. 

The ensembles of gluon field configurations were gen- 
erated by the MILC collaboration p[ 0. Different en- 
sembles include the effect of zero, two and three flavors 
of sea quarks. Ensembles with zero flavors of sea quarks 
are called 'quenched', and those with sea quarks included 
are called 'unquenched' (or 'dynamical'). The ensem- 
bles with three flavors arc most interesting from a phys- 
ical point of view because sea u, d, and s quarks are 
an important component of the vacuum of QCD in the 
real world. The MILC collaboration has made ensembles 
with a range of masses for the sea u and d quark masses 
(which are taken to be equal) down to much lighter val- 
ues than before and within reach of their real world val- 
ues. The ensembles in which the sea u and d quarks are 
much lighter than the sea s quark will be called '2+1 
flavor' ensembles. The inclusion of light sea quarks has 
been made possible by a formulation for light quarks on 
the lattice called improved staggered (asqtad) quarks 0. 
Apart from being numerically very fast, this formulation 
of the QCD action for light quarks is also very accu- 
rate. It has all sources of discretisation error removed 
to ©(a^a^), where a is the lattice spacing and as is the 
QCD coupling constant. The gluonic piece of the ac- 
tion is also simulated with a very accurate discretisation 
of QCD which has errors at 0(0^0^) apart from effects 
from sea quarks, expected to be small, which feed in at 
0{asa^) In order to check the effect of discreti- 

sation errors, two sets of ensembles are available with 



The sets of ensembles of gluon field configurations on 
which we have calculated bottomonium correlators are 
given in Tabled The ensembles within the 'coarse' set 
and the 'fine' set are matched to have approximately 
the same lattice spacing to separate discretisation effects 
from effects of changing the sea quark mass. There are 
results for a variety ofu/d quark masses, heavier than the 
real world values and we will examine the dependence of 
our results on the u/d quark mass, aiming for conclusions 
that are relevant to the physical (chiral) limit where the 
u/d quark mass is very small. The sea s quark mass 
takes only one value here, which is only approximately 
tuned (for example, it is slightly high on the coarse set 
of ensembles). Again, whether this has any impact or 
not depends on how sensitive results are to the sea quark 
mass. The b quarks are treated only as valence quarks 
because their effect in the sea should be negligible, b eing 
suppressed by inverse powers of the b quark mass [ll|. 
We also neglect c quark vacuum polarization effects for 
the same reason. 



For the valence b quarks in the T system we use lat- 
tice Nonrelativistic QCD (NRQCD) which has been de- 
veloped over many years [l3. Il3l Il4| to handle well the 
physics of heavy quark systems on the lattice. The main 
point here is that the b quark has a mass which is larger 
than 1 in lattice units and therefore momentum scales 
of order the mass cannot be simulated without large dis- 
cretisation errors. On the other hand such large scales do 
not need to be simulated on the lattice because they are 
irrelevant to the internal dynamics of the bound state 
which sets the mass splittings. The b quark is non- 
relativistic inside its bound states (w^ 0.1 for the T) 
and so a non-relativistic expansion of the QCD action is 
appropriate that accurately handles scales of the order of 
typical momenta and kinetic energies inside these states. 
Non-relativistic QCD can be matched to full QCD order 
by order in the expansion in and Ofg. 



The lattice NRQCD Hamiltonian that we use is given 
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TABLE L Parameters and details of MILC configurations used for T correlator calculations, is the number of sea flavors, 
varying from (quenched) through 2 (degenerate u and d only) and 3 (degenerate u, d and s) to the most physical 2+1 
(degenerate u and d and more massive s). 13 is the bare gauge coupling parameter in the gluon action, IQ/g^ . amu arua are 
the light sea u/d and s quark masses respectively in lattice units (using the MILC convention for quark mass, see [J), uql is 
our estimate of the maximal trace link in lattice Landau gauge used to tadpole-improve fields in our NRQCD action. aM^ is 
the bare b quark mass in lattice units, riconf is the number of configurations from the MILC ensemble used in our calculation 
with riorig origius for the heavy quark propagators per configuration. *For 92 of the 159 configurations, only 12 origins were 
used instead of 16. Tables IVlllI IIXI and Pcl give our lattice spacing determinations for these ensembles; they are approximately 
0.17fm for the super-coarse ensemble, 0.12fm for the coarse ensembles and O.OOfm for the fine ensembles. The pion masses 
corresponding to the u/d quark mass on the most chiral (i.e. lightest u/d mass) fully unquenched ensemble in each set are 270 
MeV on the super-coarse 360 MeV on the coarse and 330 MeV on the fine ensembles 0- 
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Here A is the symmetric lattice derivative and A'**) is 
the lattice discretisation of the contimmm J^i^f- ^'^t 
is the bare b quark mass. The spin-independent terms 
contribute to the (spin-averaged) differences in mass be- 
tween radial and orbital excitations in the spectrum and 
the ground-state. Spin-independent terms arc included 
through next-to-leading order above with the Darwin and 
relativistic corrections. We wiU now discuss remain- 
ing sources of systematic error in turn in the radial and 
orbital splittings and estimate their size. 



A. Systematic errors 

Relativistic errors. Errors from missing higher orders 
in the relativistic expansion should be 0(w"') = 1% m 
the radiative and orbital splittings in principle. However, 



this docs not allow for the fact that the error is set by 
the difference in the expectation value of the relativistic 
correction in the two states that make up the splitting. 
Estimates of the expectation value of powers of p in a 
potential model show that the 25* and IS" states are sim- 
ilar enough, and have similar enough expectation values, 
that the error in the difference between these states is sig- 
nificantly less than the naive expectation. We take the 
systematic error from relativistic corrections to be 0.5% 
for the 25* - 15 sphtting and 1% for the the IP - 
splitting here. These sources of error are tabulated in 
Table in 

Radiative errors. We must also consider missing radia- 
tive corrections to the terms that we are including in the 
NRQCD action. These are represented by the coefficients 
Ci that are required to match full QCD by taking ac- 
count of gluon radiation above the lattice cut-off, missing 
from the lattice theory. The coefficients take the form 



1- 



, where is a function of the bare heavy 
quark mass in lattice units, M^a. Calculation of the c!;p 
is in progress for this action. They were calculated earlier 
for a similar but different NRQCD (and gluon) action 
and found to have a value around 0.5 for aM° > 1. Here 
we have set the Ci to 1, giving expected systematic errors 
of 0{v'^as) = 2 — 3% in spin-independent splittings from 
a very naive analysis. A more sophisticated analysis, as 
for the relativistic corrections, uses the expectation value 
of in the different states from a potential model and 
calculates 2as(5(j3*)/8M^, where 5{p'^) is the difference 
in values appropriate to the splitting and the factor 
of 2 comes from having 2 h quarks in an T. The error 
estimates that this gives are shown in Table IhI The ef- 
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TABLE II: An estimate of systematic errors in the 25 — 15 and IP — 15 splittings in the T in our lattice QCD calculation 
arising from missing higher order relativistic and radiative corrections to the NRQCD action that we use. 
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TABLE III: An estimate of systematic errors in the 25 — 15 and IP — 15 splittings in the T in our lattice QCD calculation 
arising from discretisation errors in the NRQCD and gluon actions. 



feet of the Darwin term. A • E, can be estimated in a 
potential model by adding a delta function of appropri- 
ate strength (V^^) at the origin. The systematic error 
from missing radiative corrections to the Darwin term 
is then 47rasV^(0)/(3M^) where V'^(O) is the square of 
the 'wavefunction at the origin'. We estimate this using 
our results in section IVIIII and tabulate the estimates in 
Table UTI We can then add all our estimates together in 
quadrature to obtain a total estimated systematic error 
from relativistic and radiative corrections and this is also 
given in Table ^ The maximum error is 2.5% for the 
IP — 15 splitting on the super-coarse lattices. 

Note that we have tadpole- improved all the gauge fields 
appearing above by dividing the gluon fields, J7^, by a pa- 
rameter uq which cancels, in a mean-field way, the effect 
of a disparity between the lattice and continuum gluon 
fields induced by the fact that the lattice field is expo- 
nentially related to the continuum field Extra tad- 
pole diagrams appear in the quark-gluon coupling on the 
lattice and these take a fairly universal form, allowing 
a cancellation by uq. uq should be a measure of the 
difference between C/^ and the unit matrix, but is not 
uniquely specified. Here (as given in Table we use 



uql, an estimate of the trace of the average field in 
lattice Landau gauge (which maximises this trace) . Sub- 
sequent determinations of the trace of the Landau link Q 
have shown that the values used on the fine ensemble are 
underestimates by about 1%. In Figure ^ however, we 
show that using a rather different value of uq based on 
the plaquette (itop = \/plaq = 0.8677 on the coarse 
2-1-1 flavor 0.01/0.05 ensemble used for comparison) has 
no discernible effect (less than our 1.5% statistical error) 
on spin-independent, radial and orbital, splittings [l7l |. 
This also lends support to the idea that the are not 
very different from 1 because they must perturbatively 
correct for the different uq factors. Without tadpole- 
improvement the Ci would be very different from 1. 

Discretisation errors. The Hamiltonian above, Equa- 
tion^ also includes terms which correct for discretisation 
errors in the leading order spatial derivative (squared) 
and temporal derivative. These are the terms multiplied 
by C5 and cg respectively. Note that both terms have 
large numerical factors in the denominator. As above, C5 
and cg are set to 1 here, and the systematic error comes 
from neglected radiative corrections. The (D{as) pieces of 
C5 and cg were calculated in [l5l | . Indeed ci and cg are not 
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independent since these terms can be combined. Both 
Cg^'' and Cg^^ were found to be 0.5, i.e. 0(1), as expected. 
The expected systematic error in spin-independent split- 
tings is then 0{as{aM°)'^v'^ /U) = 5% x (a^inGeV-^) 
and 0{as{aMlj)v'^/8n) < 1% x (ainGeV) from a naive 
analysis. The more sophisticated analysis as above uses 
estimates of in a potential model and the results are 
tabulated in Table ITTll Corrections which arc higher or- 
der in a will have an even smaller effect because the addi- 
tional powers of aMj^ will come with additional powers of 
v^. The Darwin term also contributes to the discretisa- 
tion errors and an estimate of this is in progress. Prelimi- 
nary results 01 suggest that it has a negligible coefhcient 
and so we do not include an error for it here. There arc, 
however, additional discretisation errors coming from the 
gluon action. These were previously estimated for unim- 
proved glue in earlier work . Because the error is pro- 
portional to the difference in the square of the wave func- 
tion at the origin, again we find a bigger systematic error 
in the IP - IS sphtting than in the 2S - IS. This cor- 
rection is now reduced by a factor of for the improved 
glue configurations we are working with here. We give 
estimates for its size in Tabic Hill The total systematic 
error estimated from discretisation effects is also given 
in this Table. The error is sizeable on the super-coarse 
lattices and dominates that from missing relativistic and 
radiative corrections. On the fine lattices the estimated 
discretisation error is only 0.4% in the 2S — IS splitting, 
about the same size as the relativistic/radiative errors. 

From our error analysis it is clear that we should use 
the 25" — IS* splitting to set the lattice spacing, despite 
the fact that it is somewhat harder to obtain precise re- 
sults. Then the systematic error above is smaller than our 
statistical error on all except the super-coarse ensemble. 
Because we have results at three different values of the 
lattice spacing we can assess how good our estimates of 
systematic errors are, and we will do that below. 

It is interesting to note that the effect of sea c and b 
quarks is to modify the coefficients of the higher order 
improvement terms in the gluon action at 0{as). 
The additional coefficients are very small, 0.025as for the 
rectangle term for artic ~ I, so from the arguments above 
we can see that their effect would be entirely negligible 
in spin-independent splittings. 

E and B are the chromoelectric and chromomagnctic 
fields. They are generally defined by standard clovcrleaf 
operators (tadpole-improved) but here we use improved 
operators E and B for the E and B operators which 
appear in the spin-dependent terms. These terms are 
only included at leading order in the Hamiltonian above 
and therefore we expect systematic errors in the spin- 
splittings of 0{v'^) = 10% from higher orders in the rel- 
ativistic expansion and 0{as) = 20% from perturbative 
corrections in C3 and C4. Previous studies showed 
that errors in the cloverleaf operators could have a 
substantial effect for the fine structure splittings, such as 
the hyperfine splitting, presumably because this splitting 
is sensitive to short distances and therefore rather large 



momenta, pa, for the quark. Here we therefore attempt 
to ameliorate this effect by correcting the cloverleaf op- 
erator in these terms for the leading discretisation errors. 
This is done by the following replacement [Tsll ( u sing the 
corrected tadpole-improvement factors from |2j|): 



+ Ul{x - ajl)F^y{x ~ ajl)Uf_,{x - ajl) 

^)F^.A^) (2) 



(A^^i')] + x(-2 



where all gluon fields are understood to be tadpole- 
improved. 



B. Smearing and Fitting 

The b quark propagators (Green's functions) are calcu- 
lated on one pass through the lattice, using the evolution 
equation: 



G{x,t+l) = (1 



^)(1 _ ^Yulix) 



2n 

(l-^)"(l-^)G(-^,t) (3) 



with starting condition: 



G'(.f,0) = (t>{x). 



(4) 



4i{x) is a smearing function which, when the quark propa- 
gators are combined into meson correlators, improves the 
overlap with particular ground and excited states for a 
better signal. This will be discussed further below, n is a 
stability parameter which is chosen to tame (unphysical) 
high momentum modes of the b quark propagator which 
might otherwise cause the meson correlators to grow ex- 
ponentially with time rather than fall. It is numerically 
more convenient and safer to fit an exponentially falling 
correlator. 

The calculations presented in this paper incorporate 
a NRQCD stability parameter value of n = 2. Physical 
results should be independent of n for reasonable values 
of n and we checked that by repeating calculations on 
the coarse 2-1-1 flavor m^/fi/ms = 0.01/0.05 ensemble 
with n = 4. Figure ^ shows that no dependence on n 
was found in the resulting T energy splittings. 

During tuning runs to fix the b quark mass we obtained 
results with lower statistics for several values of A/" vary- 
ing by 0(10%). We found no discernible difference in 
IP— IS* or 25*— IS* splittings within our statistical/fitting 
errors. This is an expected feature of radial and orbital 
splittings in the heavyonium sector, since these splittings 
are observed experimentally to vary very little between 
charmonium and bottomonium. It is one of the reasons 
why these splittings arc useful to set the lattice spacing 
for an ensemble of configurations, because it can be done 
reliably without a precise tuning of the b quark mass. 
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We use 4 different smearing functions for S'-wave 
states. One is a local smearing or delta function, 
4'ioc{x) — 5xfi where is the source point. The other 
3 smearings are set up as simple functions of radial dis- 
tance r from the source point. This requires that the 
gluon fields are fixed to Coulomb gauge. Then a wave- 
function picture makes sense and we do not have to insert 
any gauge links into the smearing function. The smear- 
ings are labelled: 



hgs{r) = cxp( ) 

T 

hs{r) = (2ao - r)exp(- — ) 

Zao 

ids{r) = (27ag-18ao?- + 2r')cxp(--^) 

oao 



(5) 
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for good overlap with the ground, first excited and second 
(doubly) excited T states. We found ao=1.0 in lattice 
units to be a good value on testing the efi^ect of smear- 
ing and have used this value on all three sets of ensem- 
bles, despite the change in lattice spacing. In fact the 
(jjd smearing gives results very similar to the 0e smearing 
and needs to be improved to be more useful. For P-wave 
states we use g and e smearing functions also based on 
hydrogen-like wavefunctions, i.e. for the g smearing we 
take X, y and z times exp{—r/2ao) in each of the 3 coordi- 
nate directions. For the e smearing for P-waves we take 
r dependence of the form (Gag — r)ea:;p(— 7'/3ao), again 
multiplied by x, y or z for each direction. We find that 
it is better to take oq = 0.5 for the P-wavc states. For 
D-wave states for the g smearing function we have used 
the same r dependence and oq value as for the P-wave 
states, but multiplied by xy, yz or zx. Instead of an e 
type smearing we simply used sources made of an appro- 
priate set of delta functions, 1 lattice spacing from the 
origin. 

We use fits to a matrix of correlators with different 
smearings at source and sink of the form: 

Gmcsonirisc, Usk^t) = ^ a(risc, k)a*{nsk, k)e~^''K (6) 

k=l 

o-{n-sc/ski k) are the (real) amplitudes for state k to appear 
in the smearings used at source and sink respectively. We 
use a constrained fitting method f23 | which allows a large 
number of exponentials to be used in the fit by constrain- 
ing the way in which these exponentials can appear based 
on physical information. 

The augmentation of with a Bayesian term, 

2 2 _ 2 I 2 Cj\ 

where 

Aprior — / , ~2 \°) 

k Pt= 



FIG. 1: Dependence of T energy splittings on change in 
NRQCD stability parameter n and choice of uq: (a) n = 2, 
uoL (b) n = 4, itoL (c) n = 2, uop. The results are from the 
0.01/0.05 2+1 flavor coarse ensemble. 2S — IS etc indicates 
the splitting between radial excitations of the T and IP — IS 
the splitting between the appropriate ^Pi state and the T. 

allows the data over the entire range of t to be fitted 
to any number of exponentials in Equation |S1 The fit 
parameters, pk, were taken as the amplitudes An, the 
log of the ground state energy ln(i?o) and the logs of the 
energy differences ln(£'„+i — P„). This choice prevents 
the fit from venturing into unphysical negative energy 
space and ensures the correct ordering of states. Initial 
fits were performed with tight widths a^j. on the lowest 
energy parameters. The resulting fit parameters were 
then used as the prior pk values in the final fits where all 
widths were relaxed to CTp^. = 1. 

Figure 121 shows the dependence of the resulting ener- 
gies on the number of exponentials used in an example 
3x3 matrix fit (to correlators from the 0.01/0.05 2-1-1 
flavor ensemble) . Notice that the fit results become inde- 
pendent of Uexp, both in terms of fit error and in terms 
of goodness of fit for n^xp > 6. This is the usefulness 
of the Bayesian approach. It is very simple to choose a 
fit result because there is no trade-off between worries of 
systematic contamination from higher order exponentials 
and the fit error. Any result with n^xp > 6 will give the 
same result. In this = 7 was chosen. Notice 

that the energies of the higher radial excitations require 
a larger rigxp to converge, reasonably enough. If the only 
interest here was in the ground state energy fits from a 
smaller Ugxp value would be as good. 

III. RESULTS 

3 X 3 or 4 X 4 matrix fits were performed to the S'-wave 
correlators, and 2x2 fits performed to the P- and D-wave 
correlators. Z?-wavc correlators have only been calculated 
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FIG. 2; T energies in lattice units as a function of number 
of exponentials used in a 3 x 3 matrix fit for results on the 
0.01/0.05 coarse ensemble, x^/dof values indicate goodness 
of fit. 



0.0082/0.082 


a£;(l^Si) 


0.27755(36) 




0.770(8) 


aE{l^Pi) 


0.6674(64) 



TABLE IV: Fitted energies in lattice units for various radial 
and orbital excitations on the super-coarse ensemble. 



on the most chiral (0.01/0.05) 2+1 flavor coarse ensem- 
ble, and only for a subset of D-wave spin states. For each 
fit, the best chosen as one which save a reason- 

able x^/dof and was large enough for fitted results to 
be independent of Ugxp - The results for the fit energies 
are given in Tables IIVI Ivl and IVII for the super-coarse, 
coarse and fine lattices respectively. It is clear that the 
fitted energies are very accurate when compared to previ- 
ous results, even from the quenched approximation |20l |. 
This is partly because of the large number of configura- 
tions, their large physical volume and the large number 
of origins that we have been able to use. It is also partly 
because of the improved fitting strategy. 

In the next two sections we will discuss the spin- 
independent aspects of the spectrum, that is the differ- 
ences in mass between different radial and orbital excita- 
tions. In principle we should be using masses to calculate 
the mass splittings that are spin-averaged over the differ- 
ent spin states for each excitation. However, we cannot 
do that for the 5- wave states since no ^5*0 state has been 
seen experimentally. For the P-wave states we can spin- 
average but it is easier to take simply the mass of the ^Pi 
state. In potential models, where the spin-spin coupling 
induces a delta-function potential, the difference between 
the ^Pi mass and that of the spin-average of the '^Po,i.2 
is expected to be zero. We will examine the limits to that 



assumption in scction lVIII In the meantime we will work 
with the ^Si and ^Pi results from the Tables. 

For completeness we also give the results for different 
values of stability parameter n and tadpole-improvement 
factor uq for the 0.01/0.05 coarse ensemble in Table IVTll 
It is clear from this table that, although changing uq does 
not affect spin-independent splittings, it does provide an 
overall shift to the energies. 



IV. DETERMINING THE LATTICE SCALE 

The lattice spacing is fixed by the choice of f3 but is 
not determined until after the simulation when a physi- 
cal quantity is matched to experiment. Once the lattice 
spacing has been set in this way it can be used to give 
dimensionful predictions for other observables. 

There should be only one value for the lattice spacing 
and therefore it should not matter which hadron mass or 
mass splitting is used in fixing it. However, it has been a 
long known problem that in fact this is not the case in lat- 
tice simulations with unrealistic (or completely missing) 
light-quark vacuum polarization. This is because differ- 
ent physical quantities probe different energy scales and 
the strong coupling constant will not run correctly be- 
tween these scales unless the sea quark content (i.e. the 
quark vacuum polarization contribution) is sufficiently 
correct. This gives ambiguities in calculating physical 
observables as it is impossible to say what the 'right' lat- 
tice spacing is. 

Determinations of the lattice spacing from different 
splittings can then be used as a guide to judge how phys- 
ical is the sea quark content of the configurations being 
used. Here, the orbital IP- IS' (i.e. l^Pi - l^Si) and ra- 
dial 25"- IS* (i.e. 2^5*1 - l^Si) energy splittings indepen- 
dently fix the lattice spacing [23 ■ The results are given 



in Tables IVIIII IIXI and |X| for the super-coarse, coarse 
and fine lattices respectively. There we list the statis- 
tical/fitting uncertainties in and also give values for 
the expected systematic errors as estimated in Sectionim 
The systematic errors will affect all the estimates of 
in the same way, so it is only the differences between 
these errors that affect any discussion of differences be- 
tween determinations. We have not included an error 
from our inability to spin-average the S-wave states. The 
error from that is set by one quarter of the error in the 
spin-splitting (since we are using the vector state instead 
of the spin-average) and the spin-splitting is already only 
10% of a typical spin-independent splitting. A 20% error 
in the spin-splitting from missing relativistic corrections 
would then lead to a 0.5% error in the determination of 
from, for example, the IP — 15* splitting. This is 
neghgible. 

It is quite clear from these tables that there is a prob- 
lem of internal consistency in the quenched approxima- 
tion. The ratio of from IP — IS to that from 
25" — IS is 1.09(2)(4) on the coarse quenched ensem- 
ble and 1.09(2)(1) on the fine quenched ensemble. The 
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0.01/0.05 0.02/0.05 0.03/0.05 0.05/0.05 



0.02 



Quenched 



aEil'^Si) 0.30506(27) 0.30118(25) 0.29848(26) 0.29588(34) 0.29563(26) 0.26621(24) 



a_E(2^5'i) 

o£(l'Pi) 

aE{V^D2) 
aE{2^D2) 



0.6578(56) 
0.887(20) 
0.5852(37) 
0.822(17) 
0.751(12) 
1.026(34) 



0.6519(44) 
0.878(22) 
0.5705(41) 
0.801(16) 



0.627(10) 
0.870(32) 
0.5593(70) 
0.784(53) 



0.6381(95) 
0.927(21) 
0.5593(61) 
0.827(18) 



0.6326(79) 
0.851(33) 
0.5519(59) 
0.787(29) 



0.6073(49) 
0.887(12) 
0.5111(27) 
0.788(12) 



TABLE V: Fitted energies in lattice units for different radial and orbital excitations for the coarse lattice ensembles. 



0.0062/0.031 0.0124/0.031 Quenched 



0.0082/0.082 



aE{l^Si) 
aE{2^Si) 
aE{3^Si) 
aE{l^Pi) 
aE(2^Pi) 



0.23248(46) 
0.4818(33) 

0.619(14) 
0.4233(37) 

0.574(44) 



0.28677(40) 
0.5303(27) 

0.671(13) 
0.4705(30) 

0.628(16) 



0.18517(23) 
0.42734(25) 
0.6472(71) 
0.3595(19) 
0.5470(86) 



TABLE VL Fitted energies in lattice units for various radial 
and orbital excitations on the fine lattice ensembles. 



,s(GeV) 1.144(19)(11)(23) 
,g(GeV) 1.128(19)(28)(90) 

TABLE VIIL Inverse lattice spacing determinations for the 
super-coarse ensembles from both the l^Pi — l^Si and 2^5*1 — 
I'^Si splittings. The first error given is statistical/fitting. 
The second is the expected systematic error from relativis- 
tic/radiative corrections and the third from discretisation er- 



plaq 



aE{l^Si) 0.29433(27) 0.47736(25) 

aE{2^Si) 0.6449(65) 0.8296(62) 

a£;(3='5'i) 0.8799(65) 1.0652(62) 

aS(l^Pi) 0.5743(37) 0.7599(45) 
aE{2^Pi) 0.811(16) 0.995(22) 

TABLE VII: T spin independent energy splittings in lattice 
units on the coarse 2+1 flavor 0.01/0.05 ensemble. These 
were used for checking systematic errors by comparison with 
results in Table [V] 



first error here is from statistics/fitting and the second 
from an estimate of the difference in systematic errors 
for the IP — IS and 25* — 15 spUttings. Previous re- 
sults from quenched NRQCD calculations also gave ra- 
tios larger than 1 |23, U3> HE lUli ahhough generally 
with lower significance because of poorer statistics or fit- 
ting. One problem with previous results was that most 
of them used an unimproved gluon action and the effect 
of the tree-level 0{a^) discretisation errors is to reduce 
this ratio of determinations by a few percent, ob- 
scuring the disagreement with 1. The ratio including a 
correction for this is given in reference j23| as 1.06(4) at 
P = 6.0. Without correcting for it, reference [2^1 gives 
1.12(9). It is clear that our new quenched results on im- 
proved glue are a considerable improvement in clarity on 
these previous results. 

Figure |31 shows the ratio of lattice spacing determi- 
nations from the 2S — IS and IP — IS splittings as a 
function of the bare light (u/d) quark mass in our sim- 
ulations (see Table . Also included are results with 
2 flavors of sea quarks only, and over to the right of 



1.15 
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0.95 



0.02 0.04 0.06 0.08 0.1 oo 

miig,,(GeV) 

FIG. 3: The ratio of lattice spacing values obtained from the 
25—15' and IP — IS splittings in the T system as a function of 
bare light quark mass in the lattice QCD calculation. Notice 
the very magnified j/-axis scale. The crossed square is from 
the super-coarse ensemble, the filled squares from the coarse 
2-f 1 flavor ensembles and the open squares from the fine 2-1-1 
flavor ensembles. The flUed diamond is from the coarse Uf — 2 
ensemble. The open and filled triangles on the right of the 
plot are results from the quenched fine and coarse ensembles 
respectively. Errors are statistical only. 



the plot, the quenched results. The ratio of inverse lat- 
tice spacings on the most chiral 2+1 flavor ensemble at 
each lattice spacing gives 0.981 (23) (60) for super-coarse, 
0.984(20)(34) for coarse and 1.021(23)(13) for fine, where 
the first error is statistical and the second systematic. All 
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0.01/0.05 0.02/0.05 0.03/0.05 0.05/0.05 002 Quenched 

a^s-isiGeV) 1.596(25)(11)(16) 1.605(20)(11)(16) 1.714(52)(12)(17) 1.645(46)(12)(16) 1.671(39)(12)(17) 1.651(24)(12)(17) 
«ip-is(GeV) 1.571(21)(31)(63) 1.634(25) (33) (65) 1.687(45) (34) (67) 1.670(39) (33) (67) 1.717(40)(34)(68) 1.797(20)(36)(72) 



TABLE IX: Inverse lattice spacing determinations for the coarse ensembles from both the l^Pi — I'' Si and 2^5*1 — l^Si splittings. 
The first error given is statistical/fitting. The second is the expected systematic error from relativistic/radiative corrections 
and the third from discretisation errors. 



0.0062/0.031 0.0124/0.031 Quenched 
a2S-is(GeV) 2.258(30) (14) (9) 2.312(26)(14)(9) 2.324(24) (14) (9) 
a~p„ig(GeV) 2.305(45)(35)(35) 2.390(38)(36)(36) 2.524(28) (38) (38) 



TABLE X: Inverse lattice spacing determinations for the fine ensembles from both the l^Pi — V^Si and 2'^5i — l^Si splittings. 
The first error given is statistical/fitting. The second is the expected systematic error from relativistic/radiative corrections 
and the third from discretisation errors. 



of these results are in agreement with 1, unlike the case in 
the quenched approximation. The results agree well with 
each other too, and this provides a check on our estimates 
of systematic errors from discretisation effects. For the 
fine lattices discretisation effects are small, but on the 
super-coarse lattices we have estimated that the IP — IS* 
spHtting could be moved by 6% more than the 25* — IS* 
splitting by discretisation and relativistic/radiative sys- 
tematic errors. This looks like an overestimate since our 
super-coarse result agrees with 1 within statistical errors 
of 2%. 

The agreement between the inverse lattice spacing re- 
sults from the two different splittings is confirmation of 
the approach to the physical real world of a single inverse 
lattice spacing as the sea quark content is made realistic. 
Some signs of this were seen in earlier 2 flavor simula- 
tions [ij, m El l23- Indeed here our 2 flavor results 
also give agreement with 1 for the ratio of lattice spacing 
values within the statistical error. As discussed earlier, 
the lattice spacing from the 2S ~ IS splitting is to be 
preferred on the grounds of systematic error and so we 
will use that in what follows. 

We can also ask what the lattice spacing derived from 
T splittings implies for other quantities. A quantity that 
is used frequently to set the scale is the ro parameter . 
It is derived from the heavy quark potential as the value 
of r at which r^F(r) = 1.65, where F{r) is the gradi- 
ent of the potential, rg is easily and precisely calcula- 
ble in lattice calculations but has the disadvantage that 
its value in physical units is not directly given by ex- 
periment. Instead its value must be given by numerical 
simulation on the lattice and comparison to an experi- 
mentally accessible quantity. Here we can give a value 
for To from comparison to T splittings. In Figure 0] we 
plot To in fm as a function of light sea quark mass, fixing 
the lattice spacing from the 25" — IS" splitting. The value 
of To in lattice units has been calculated by the MILC 
collaboration H Q- The result does not show any de- 
pendence on either the light quark mass for light enough 
light quark mass or the lattice spacing. Our best result 



is from the fine lattices, and gives a physical value for tq 
equal to 0.469(7)fm. The error includes the systematic 
errors in the determination of the lattice spacing on the 
fine lattices and a similar systematic error for the deter- 
mination of tq . However these errors are small compared 
to the statistical error from the determination of the lat- 
tice spacing. This result is in agreement with the results 
from which were based on our preliminary analysis 
of the T results presented in this paper. We similarly 
have a value for the heavy quark potential variable ri 
(where r^F{r) = 1.00) of 0.321 (5)fm. This is also 
shown as a function of light quark mass in Figure [S] and 
also agrees with the preliminary value from 01 . This vari- 
able in fact shows a very flat curve as a function of light 
quark mass when the scale is set by the T 25* — IS" split- 
ting, i.e. the light quark vacuum polarisation effects are 
very similar in both variables. Note that if we had used 
the IP — IS splitting to set the lattice spacing, both ri 
and To would have had values in fm differing by 10% in 
the quenched approximation, because of the ambiguity 
in setting the lattice spacing demonstrated in Figure 13 
The result on the chiral 2-1-1 flavor configurations would 
be unchanged, because here there is no ambiguity about 
the value of a 0. 



V. COMPARING RADIAL AND ORBITAL 
SPLITTINGS TO EXPERIMENT 

The 25* — IS" inverse lattice spacing values and fit ener- 
gies from the previous sections were combined to give di- 
mensionful determinations of the radial and orbital split- 
tings 3^Si - l^Si, l^Pi - l^Si, 2iPi - l^Si. FigureEl 
plots the results as a function of the sea light quark mass 
for the coarse and fine ensembles. Bursts denote the ex- 
perimental values where we have used the spin-average 
of the Xb states in place of the unobserved h(,. 

Good agreement with experiment is seen as the chi- 
ral limit is approached and the flatness of the approach 
indicates that chiral extrapolations would have no signif- 
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>"light(GeV) 

FIG. 4: The value of ro in fm obtained on diflerent super- 
coarse, coarse and fine ensembles, using the 2S — IS splitting 
in the T system to fix the lattice spacing. The crossed square 
is from the super-coarse ensemble, the filled squares from the 
coarse 2-1-1 flavor ensembles and the open squares from the 
fine 2-1-1 flavor ensembles. The open and filled triangles on the 
right of the plot are results from the quenched fine and coarse 
ensembles respectively. The j/-axis scale is expanded to make 
error bars visible, ro values in lattice units are taken from 0, 
0, ■ Errors include the statistical errors from determining 
the lattice spacing and from determining rg. 




0.02 0.04 0.06 0.08 0.1 
>"light(GeV) 

FIG. 5: The value of ri in fm obtained on different super- 
coarse, coarse and fine ensembles, using the 2S — IS splitting 
in the T system to fix the lattice spacing. The crossed square 
is from the super-coarse ensemble, the filled squares from the 
coarse 2-1-1 flavor ensembles and the open squares from the 
fine 2-1-1 flavor ensembles. The open and filled triangles on the 
right of the plot are results from the quenched fine and coarse 
ensembles respectively. The j/-axis scale is expanded to make 
error bars visible, ri values in lattice units are taken from 0, 
0, 13 ■ Errors include the statistical errors from determining 
the lattice spacing and from determining ri. 



icant effect in shifting the results from the most chiral 
ensemble. 

Figure [7| shows the radial and orbital levels in the T 
spectrum with the lattice results from the coarse and fine 
2-1-1 flavor ensembles with lightest sea quark mass con- 
trasted to coarse and fine quenched results. The solid 
lines are the experimental values, including the recently 
observed state 0. On the lattice the spin 2 rep- 
resentation is split into two representations labelled E 
and T. Here we used the E representation. This will 
be discussed further in section IV 111 The 15 energy and 
2S ~ IS splitting arc not predictions as they have already 
been used to fix the b quark mass in lattice units and the 
lattice spacing respectively. The quenched results clearly 
do not agree with experiment for other splittings whereas 
the 2-1-1 flavor results do. 



VI. THE b QUARK MASS 

The bare b quark mass in the lattice Lagrangian is 
fixed by the requirement to get the T mass correct. Be- 
cause the quark mass term has been removed from the 
NRQCD Lagrangian the energies calculated for hadrons 
at zero momentum are shifted from their masses. This 
means that energy differences (for particles containing 
the same number of quarks or antiquarks) do give mass 
differences but in order to fix one mass absolutely, a dis- 



persion relation must be calculated. In sections HVl and IVI 
we discussed energy differences and mass splittings. Here 
we give results for the mass of the T from calculations of 
its energy as a function of lattice momentum, and discuss 
how well the b quark mass can then be determined. 



A. Fixing the T mass 

The continuum relationship between energy, E, and 
spatial momentum, p, for a hadron in a theory such as 
NRQCD where the zero of energy has been shifted is 
given by: 

E{p) = £:(0) + Vp^TaP _ M. (9) 

The hadron mass is denoted by M here, sometimes called 
the kinetic mass to emphasise that it is determined from 
a dispersion relation rather than E{0). Equation is a 
fully relativistic formula. When used in our lattice calcu- 
lations, p is a lattice momentum given in lattice units by 
2'K{nx,ny,nz) / L where L is the number of lattice spac- 
ings in the spatial directions (16 for super-coarse, 20 for 
coarse and 28 for fine) and ?^z:,y,z take integer values be- 
tween and L — 1. Mkm is obtained by first fitting the 
energy of jack-knifed ratios of T correlators at zero and 
non-zero values of p to an exponential form. The expo- 
nent yields the energy difference AE = E{p) — E{0) and 
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FIG. 6: Spin-independent splittings in the T spectrum as a 
function of bare sea light quark mass. The IP — IS, 2P — \S, 
and 35 — 15 splittings are shown. The S state is always the 
^5*1. The P state from the lattice calculations is the ^P\ \ from 
experiment it is the spin-average of ^Po,i,2- Closed squares 
are from coarse 2-1-1 flavor lattices and open squares from 
fine 2-f 1 flavor lattices. The crossed square for the IP — IS 
splitting is from the super-coarse ensemble. Closed diamonds 
are from the coarse n/ = 2 run. Closed and open triangles are 
from quenched coarse and fine ensembles respectively. Bursts 
give the experimental results. Errors include the statistical 
errors from determining the lattice spacing from the 2S — \S 
splitting. 



M 



2AE 



(10) 



Note that a non-rclativistic formula would not include the 
(A£')^ term, and this becomes increasingly negligible as 
the mass increases. 

Table Ixll gives . as an example of our results, values for 
AE and M obtained on the coarse 2 -f 1 flavor 0.01/0.05 
ensemble for a range of lattice momenta. These results 
are plotted in physical units (using the lattice spacing 
from the 25*— IS* splitting of section llV)l in Figure|Hl Fig- 
urelHlshows that the T mass is stable and well-determined 
out to very large values of p^, far larger than are rou- 
tinely accessible in light hadron calculations. The value 
does not change significantly, showing that lattice arte- 
fact terms and errors coming from missing relativistic 
corrections to the NRQCD action in the lattice disper- 
sion relation are not significant. 

An equivalent illustration of the stability of the T mass 
is given by the ratio of the p and E terms in Equation |51 
This ratio is called the 'speed of light' and should have 
value 1 in our units. We calculate it as: 



(AE)^ + 2MAE 



(11) 
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FIG. 7: The T spectrum of radial and orbital levels, taking the 
most chiral points from Figure El Closed and open symbols 
are from coarse and fine lattices respectively. Squares and tri- 
angles denote unquenched and quenched results respectively. 
Lines represent experiment. 





aAE 


aA/x 




1 


0.008264(40) 


5.967(29) 




2 


0.016507(84) 


5.971(30) 


0.999(5) 


3 


0.02473(13) 


5.974(32) 


0.999(5) 


4 


0.03286(19) 


5.990(34) 


0.996(6) 


5 


0.04104(25) 


5.991(36) 


0.996(6) 


6 


0.04928(31) 


5.983(37) 


0.997(6) 


8 


0.06540(47) 


6.004(43) 


0.994(7) 


9 


0.07353(57) 


6.003(46) 


0.994(8) 


12 


0.09778(92) 


6.007(56) 


0.993(9) 



TABLE XI: Energy differences of zero and non-zero momen- 
tum T mesons, and resulting T masses, in lattice units. Non- 
zero momenta are denoted by = ri^ + Uy + n1. These are 
unambiguous except for = 9 which corresponds to the n- 
vector (2,2,1) and its permutations. The fourth column gives 
values for the speed of light, c^, defined in the text. The 
results are for the coarse 2-1-1 fiavor 0.01/0.05 ensemble. 



where Af is determined from above using values of p cor- 
responding to rt^ = and 1 . This is a quantity much used 
as a check of how well discretisation errors are controlled 
in actions for light quarks on the lattice [s^l- Even the 
best light quark actions show discrepancies of from 1 of 
a few percent at typical lattice spacing values; some of the 
poorer actions show discrepancies of many percent |3ll |. 
Our results are tabulated in Table Ixll and plotted in Fig- 
ure 1^ They show that never deviates from 1 by as 
much as 1% and we can determine it with errors of less 
than 1%. 

Deviations from 1 of would be caused by both miss- 
ing higher order relativistic corrections in the NRQCD 
action and discretisation errors in the NRQCD and gluon 
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FIG. 8: The value of the T mass in physical units plotted 
against the the lattice momentum for the n/ — 2 + 1 0.01/0.05 
coarse ensemble for aM^ — 2.8. Errors included are only those 
from determining Mr, since the errors from determining the 
lattice spacing are the same for all points. 



FIG. 9: The 'square of the speed of light' (see text for def- 
inition) from T correlators at non-zero momentum for the 
Uf ^ 2 + 1 0.01/0.05 coarse ensemble and aM° = 2.8. Note 
the large magnification of the y-axis scale to make the errors 
visible. 



actions. The terms that give rise to these errors have 
been discussed in Sectional Here we can estimate their 
effect on by assuming that the momentum of the T 
is shared between the two quarks equally. The rela- 
tive error in AE from missing radiative corrections to 
the (A^^))^ relativistic correction in Equation is then 
as(p/2)2/(2A4°)2. This only amoimts to 0.5% even at 
momentum (2, 2, 2) on the coarse set of ensembles, and 
0.4% on the fine set. Similarly missing radiative correc- 
tions to the discretisation errors, i.e. the effect of setting 
C5 and cg to 1, amount to no more than 0.5% at mo- 
mentum (2,2,2) on the coarse set of ensembles, 1% on 
the super-coarse set and 0.2% on the fine set. These 
estimates bear out the results shown in Figure El and em- 
phasise how well the NRQCD formalism handles heavy 
quarks at low non-zero momenta. This is important not 
only for extracting hadron masses but also for the analy- 
sis of form factors for B meson semi-leptonic decay |32l |. 

We obtain values for the T mass from Equation [TUI 
using momenta and = 1. These are given in Ta- 
bles l^n] for our super-coarse, coarse and fine 
ensembles respectively. 



B. Determining the b quark mass 

The bare b mass. The values for the T masses are in 
agreement with the experimental T mass for most of our 
ensembles, indicating that the bare quark masses used for 
the b quarks are good ones. There is particularly close 
agreement for the most important ensembles - the most 
chiral coarse and fine unquenched ensembles. Because 
the calculated T mass in physical units is proportional to 



0.0082/0.082 



aMr 8.086(17) 
Mr /GeV 9.25(15)(21) 



TABLE XII: The T mass in lattice units and physical units 
obtained from momenta and — 1 for the super-coarse en- 
semble. The results in lattice units include statistical/fitting 
errors only; the results in physical units include two errors. 
The first is statistical/fitting, dominated by that from deter- 
mining the lattice spacing from the 25 — IS splitting. The 
second error is the combined systematic error from the lattice 
spacing determination. 



a good approximation to the bare quark mass in physical 
units, shifts of the bare lattice quark mass and shifts 
of the lattice spacing compensate each other [3^ when 
calculating the bare quark mass in physical units. This 
means that errors in the lattice spacing determination 
or in the tuning of the bare quark mass in lattice units 
do not in fact feed through into errors in the bare quark 
mass in GeV. This is demonstrated most clearly if we 
determine 



M- 



T ,expt 



'aMr 



latt 



(12) 



where aM? . is the bare lattice b quark mass in lattice 

o Sim ^ 

units used in the calculation (see Table HI). Wc tabulate 
the values of M^{a) from this formula for each of our lat- 
tice spacings in Table IXVl The errors quoted include sta- 
tistical errors (which range from 0.2% to 0.5%) as well as 
systematic errors from relativistic/discretization errors. 
The latter are estimated by noting that the binding en- 
ergy is only 6% of the T mass (as we shall see below). 
The leading errors, from missing radiative corrections to 
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0.01/0.05 


0.02/0.05 


0.03/0.05 0.05/0.05 


0.02 


Quenched 


aMr 5.968(29) 
Mr/GeV 9.52(16)(11) 


5.942(29) 
9.54(13)(11) 


5.935(26) 5.906(34) 
10.17(31)(12) 9.72(28)(11) 


5.875(27) 
9.82(23)(11) 


5.855(37) 
9.67(15)(11) 



TABLE XIII: The T mass in lattice units and physical units obtained from momenta and = 1 for the coarse set of 
ensembles. The results in lattice units include statistical/fitting errors only; the results in physical units include two errors. 
The first is statistical/fitting, dominated by that from determining the lattice spacing from the 25* — 15' splitting. The second 
error is the combined systematic error from the lattice spacing determination. 



0.0062/0.031 0.0124/0.031 Quenched 

aMr 4.218(18) 4.268(15) 4.145(12) 
Mr/GeV 9.52(13)(7) 9.87(12)(7) 9.63(10)(7) 



TABLE XIV: The T mass in lattice units and physical units 
obtained from momenta and — 1 for the fine set of en- 
sembles. The results in lattice units include statistical/fitting 
errors only; the results in physical units include two errors. 
The first is statistical/fitting, dominated by that from deter- 
mining the lattice spacing from the 2S — IS splitting. The 
second error is the combined systematic error from the lattice 
spacing determination. 

the leading relativistic/discretization corrections, are ex- 
pected to be 20-30% of that 6%. Note that the bare 
mass, M^(a), decreases with decreasing lattice spacing, 
as expected for a running mass. The masses from the 
comparable quenched simulations are only slightly higher 
at 4.52 GeV on the coarse ensemble and 4.45 GeV on the 
fine ensemble. We now discuss the conversion of the bare 
mass first to a perturbative pole mass and then to an 
MS mass for the b quark. 

Renormalization of the b quark mass. There are two 
methods for determining the b pole mass, defined per- 
turbatively, from an NRQCD simulation [33. Both use 
perturbative results obtained from the heavy quark self- 
energy [s^l- The first uses perturbation theory to relate 
the pole mass directly to M°(a), defined above: 

Af = Z,n MS = (1 + Cm «. + ••• ) M° (13) 

where Cm ^depends on ail/" and values are given in Ta- 
ble IX VI (34i . The second method uses perturbation the- 
ory to relate the NRQCD energy of the T to the 'binding 
energy' of the meson and therefore to the pole mass: 

Mr'" - liMr^e^pt - iE,,m - 2Eo)) (14) 

where Esi^ is the fitted NRQCD energy from T correla- 
tors at zero momentum given in Tables IIVI Ivl IVII Here 
Eq is the NRQCD energy of an isolated b quark, com- 
puted in perturbation theory: 

aEo = CEoas + ■■■ . (15) 

Ceq again depends upon aMjj and the appropriate values 
are given in Table |XV)| . 




2 4 6 8 10 



non-pert shift 

FIG. 10: The perturbative estimate of the shift between mass 
and simulation energy for the T versus the result obtained 
non-perturbatively from our calculations. The results are 
given in lattice units for the most chiral 2-1-1 flavor ensem- 
bles on super-coarse (0.0082/0.082), coarse (0.01/0.05) and 
fine (0.0062/0.031) lattices (right to left). The errors on the 
perturbative points are estimates of unknown terms, taken 
as 2 X «y(g*) for 2 b quarks in an T; the errors on the non- 
perturbative results are negligible in comparison (0.5%). The 
line drawn is y = x for comparison. 



Unfortunately the momentum scales for as in our two 
perturbative formulas are quite small. Working in the 
V scheme and using BLM scale fixing [l6l l3^ an earlier 
analysis, using a slightly different version of NRQCD and 
an unimproved gluon action, found scales of order 0.6/a 
for both the Zm and aEo expansions uM- Even on our 
finest lattices, Q;y(0.6/a) is about 0.75 [33. The small 
scale is expected. It refiects real infrared sensitivity in the 
pole mass. Such large values for the coupling constant 
mean that we cannot obtain an accurate estimate of the 
pole mass. 

By equating our two formulas for M^°'°, we obtain an 
interesting relationship between b quark and T-meson 
properties: 

2(ZmA/0 - £^0) = Mr - Ssim (16) 

The quantities on the right-hand side of this equation 
are measured nonperturbativcly in our simulations, while 
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Cm 


Ceo 




super-coarse 


4.68(9) 


0.082 


0.850 


4.31(44) 


coarse 


4.44(7) 


0.235 


0.767 


4.30(31) 


fine 


4.37(5) 


0.421 


0.689 


4.44(26) 



TABLE XV: The b quark mass in physical units from the most chiral 2+1 flavor lattice ensembles (i.e. 0.0082/0.082 super- 
coarse, 0.01/0.05 coarse and 0.0062/0.031 fine). The bare mass is given on the left, the perturbative coefficients in Zm and Eo 
(see text) in the centre and the mass in the MS scheme at its own scale on the right. For the bare mass the error is dominated 
by a systematic error from the NRQCD action; for the MS mass the error is an estimate from unknown higher orders in 
perturbation theory. 



the left-hand side is determined from perturbation the- 
ory. The infrared sensitivity that plagues the pole mass 
cancels between the two terms on the right-hand side, so 
that scales q* for ay vary from approximately 1.7/a to 
1.3/a from super-coarse to fine lattice spacingsj^. We 
compare the perturbative predictions (left-hand side) and 
our nonperturbative measurements (right-hand side) for 
this shift in Figure EH The agreement is excellent even 
as this shift varies by almost a factor of two. 

The MS b quark mass. Our inability to obtain a pole 
mass is not much of a handicap since most applications 
require the MS mass. The MS mass, like A'/°, is a bare 
mass, and therefore is free of the infrared problems asso- 
ciated with the pole mass. We can convert our previous 
formulas for M^°'® into formulas for using the stan- 
dard formula, 

^n^) = (i - V + ■ ■ ■ ) 

(17) 

Combining this formula with Eg nation 1131 gives 

Mp(/i) = (l + ^ {Cmtt - ^ - 2 log(^)) + . . . ) il4". 

_ ' (18) 

The choice of the MS scale /i is arbitrary, but we want 
to choose a scale that matches the MS and lattice cut- 
offs — so that the two bare masses, Af^^(/i) and are 
roughly equivalent. An obvious choice is the scale 
such that the first-order coefficient in Equation 1181 van- 
ishes; this gives values for a/i* between 2.4 and 1.9 for our 
three lattice spacings, showing only weak a-dependence 
(as expected when aMf, 1). We simphfy the analysis 
here by setting fi = 2/ a for all of our lattice spacings, 
resulting in first-order coefficients that range from 0.16 
to 0.06. Previous work [s^ suggests that an appropri- 
ate coupling constant in this formula is av(2/a), where, 
as expected, the scale is much larger than for the pole- 
mass formulas. We again take values for ay from |35| . 
Our final results, evolved to the MS mass using 3-loop 
evolution ; are quoted in Table IXVl These include a 
perturbative error uncertainty which we estimate to be 
1 X ay(2/a). Our results from different lattice spacings 
are in excellent agreement; we take 4.4(3) GeV as our 
final result. 

Further checks of the b quark mass. An important 
check of the b quark mass is whether the masses of other 



hadrons containing b quarks calculated on the lattice 
agree with experiment when the b quark mass is fixed 
from the T mass. In the quenched approximation this 
has always been problematic as a result of the ambiguity 
in determination of the lattice spacing. This led to large 
differences in b quark masses used in B calculations and 
T calculations - 30% at /3=6.0 [Tll40l|. In a lattice sim- 
ulation which matches the real world there is only one 
value for the lattice spacing and one value for the bare b 
quark mass. Here we test this for the MILC 2-1-1 flavor 
configurations. 

Figure [TTl gives the mass for the Bg meson on the 2-t-l 
flavor 0.01/0.05 coarse ensemble using a bare b quark 
mass that is the same as that of the T calculation and 
flxing the lattice spacing from the 25' — 15 T splitting, 
again as in the T calculation. The Bg mass is obtained 
from the dispersion relation. Equation |51 as for the T. 
The s quark mass is flxed from the K mass The 
experimental Bg mass is 5.37 GcV so the lattice results 
are 4% or Icr high, i.e. in agreement within errors. This 
is a big improvement over the results in the quenched 
approximation and shows that these unquenched config- 
urations are at last yielding an unambiguous picture of 
the real world. 

There is a potential systematic error of 1-2%, as ex- 
plained above, coming from the NRQCD action in the T 
case. As we now explain, this source of error is not mir- 
rored in the B system, and so does remain a potential sys- 
tematic difference between the two. This can certainly 
explain the results in Figure 1111 if they are not simply 
statistical fluctuation, and would indicate that further 
improvements to the NRQCD action would remove any 
discrepancy. 

The power-counting for the NRQCD action for a 
heavy-light meson is quite different from that of a heavy- 
heavy meson, since it is then powers of 1 /Mt that matter, 
with Aqco/Mb « 10% for the B. All terms through 
second-order (1/M^) are included at tree- level in the 
NRQCD action we have used. Equation^ and we expect 
smaller systematic errors than in T system. Formally 
the largest systematic error comes from radiative correc- 
tions to the (J ■ B/Mb term but this term itself does not 
have a large effect as can be seen from comparing the B 
meson hyperfine splitting to the B mass. The B meson 
binding energy {Egim — Eq) is quite a different size, and 
of different sign, to that of the T. It arises largely from 
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p (GeV) 

FIG. 11: The value of the Bs mass in physical units plotted 
against the lattice momentum for the 71/ —2 + 1 0.01/0.05 
coarse ensemble for aM° — 2.8, and anis — 0.04 using the 
MILC convention for light quark mass. Errors included are 
only those from determining Mki„, since the errors from de- 
termining the scale are the same for all points. 



the light quark sector and here errors from the improved 
staggered action appear. These are formally as (aAgcD)^ 
which can be estimated to be around 2%, but applied to 
the binding energy this gives a very small error in the B 
meson mass. 

Another check of agreement with experiment with 
unquenched lattice QCD is to study 2mB — m-f or 
'^TTiB^ — rriY- Explicit terms in the b quark mass can- 
cel in this difference leaving sensitivity to the differences 
in binding energies in the two systems (with slight heavy 
quark mass dependence). From the discussion above we 
expect systematic errors of order 2% in this quantity, 
coming either from errors in the NRQCD action on the 
heavy-heavy side or from discretisation errors in the light 
quark action on the heavy-light side. To avoid errors on 
the heavy-light side coming from radiative corrections to 
the a-B term, we spin-average over pseudoscalar and vec- 
tor heavy-light states. The results are shown in Figure [T^ 
and show indeed good agreement with the experimental 
results at the level expected. This is again strong confir- 
mation that lattice QCD with light sea quarks reproduces 
experimental hadron masses across a wide range of the 
spectrum. 
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FIG. 12: The mass difference between twice the spin-average 
of B and B* and the T for B^/a and Bs, plotted against the 
sea u/d mass in physical units for coarse unquenched config- 
urations. The results are for the 2-1-1 flavor 0.01/0.05 and 
0.02/0.05 lattices. The u/d mass used in the Bu/d is then 
0.01 or 0.02 as appropriate and the s mass used in the Bs 
is 0.04, fixed from the K mass 0. The upper points are for 
the Bs and the lower for the -B„/^. The experimental results 
are given by the bursts |2^ where we have assumed that the 
Bg-Bs splitting is the same as that for the B^^^-B^^/^. Errors 
include statistics/fitting and errors. 



0.0082/0.082 


aE{2^Si) - 

£{2^Si)- 


-aE{l'So) 
- aE{2^So) 

-B(2lSo) 


0.0348(5) 
0.018(13) 
0.52(38) 



TABLE XVI: S-wave fine structure in the T spectrum for the 
super-coarse 2-1-1 flavor lattice ensemble. Energies are given 
in lattice units. Errors are statistical/fitting only. 



obtain are only accurate to 0{v'^) ~ 10%. In addition 
there are 0{as) corrections to the coefficients of these 
terms, such as C4 in equation ^ which could amount to 
a 20-30% error. In fact, as we shall see, we are able to 
do rather better than this, but results in this section are 
still rather qualitative. We present them nevertheless as 
providing useful pointers for future calculations. 



S-wave fine structure 



VII. FINE STRUCTURE IN THE T SPECTRUM 

Fine structure in the T spectrum arises from terms in 
the action which are suppressed by a power of v'^ com- 
pared to the leading term, Hq. Such sub-leading terms 
themselves have sub-leading corrections which we have 
not included. We expect therefore that the results we 



There are no experimental results for either the ground 
or radially excited S'-wave spin splitting (hyperfine split- 
ting) in the T system, and so these are key results to be 
predicted by lattice QCD. This splitting is the mass dif- 
ference between the vector (T) state and the pseudoscalar 
{rjb) and so a prediction of the splitting is equivalent to 
a prediction of the unknown yy^''' mass. Our result for 



16 



0.01/0.05 0.02/0.05 0.03/0.05 0.05/0.05 0.02 Quenched 

aE(V^Si)-aE{l^So) 0.03128(34) 0.03062(32) 0.02972(33) 0.02925(43) 0.02668(34) 0.02152(30) 

aE{2^Si)-aE{2^So) 0.0111(75) 0.0172(62) 0.017(14) 0.013(13) 0.0102(109) 0.0137(78) 
fiqag'jllliis"! 0.35(24) 0.56(20) 0.57(47) 0.44(44) 0.38(41) 0.64(36) 

TABLE XVII: 5-wave fine structure in the T spectrum for the coarse ensembles. Energies are given in lattice units. Errors are 
statistical/fitting only. 







0.0062/0.031 0.0124/0.031 


Quenched 


aE{l^Si) 
aE{2'^Si) 

E(l^Si 


-aE(l^So) 
~aE{2^So) 

-E(2lSo) 
-E(llSo) 


0.02599(53) 0.02422(45) 
0.0141(55) 0.0120(42) 
0.54(21) 0.50(18) 


0.01735(29) 
0.0107(31) 
0.62(18) 



TABLE XVIII: S'-wave fine structure in the T spectrum for 
the fine lattice ensembles. Energies are given in lattice units. 
Errors are statistical/fitting only. 




0.02 0.04 0.06 0.08 0.1 
ni|igh,(GeV) 

FIG. 13: The ratio of the splitting between the T' and rjl 
to that between the T and the rjb. Crossed squares are from 
super-coarse lattices, closed squares are from coarse lattices 
and open squares from fine lattices. The closed diamond is 
from the coarse n/ = 2 ensemble. The closed and open tri- 
angles are from the coarse and fine quenched ensembles. The 
burst represents the current experimental result for the sys- 
tem, using the value for the mass of the rj'^ from the BELLE 
collaboration 1411. 



the ground-state hyperfine splitting is 61(14) McV, giv- 
ing M^i^ — 9.399(14) GeV. Our prediction for the radially 
excited hyperfine splitting is 30(19) MeV, giving A/,,^ = 
9.993(19) GeV. The largest contribution to the error for 
the ground-state splitting comes from unknown radiative 
corrections to the spin-dependent term in the NRQCD 
action that gives rise to this splitting. For the excited 
state splitting the dominant error is statistical. How we 
estimate the size of the errors and how we can reduce 
them in future is described in more detail below. 




0.06 0.08 0.1 
miight(GeV) 

FIG. 14: The ground-state hyperfine splitting in the T system 
as a function of the light quark mass included in the quark vac- 
uum polarization. Crossed, closed and open symbols are from 
super-coarse, coarse and fine lattices respectively. Squares 
and triangles denote unquenched and quenched results re- 
spectively. The closed diamond is from the coarse n/ = 2 
ensemble. Errors include both statistical/fitting errors from 
determining the splitting and from determining the lattice 
spacing scale. This latter error appears doubled (see text). 
Note that there is an additional overall systematic error of 
25% discussed in the text. 



Results for the ground and radially excited hyperfine 
splitting in lattice units, and their ratio are given in Ta- 
bles IXVII IX VIII and IX Villi for the super-coarse, coarse 
and fine MILC lattices respectively. Splittings are ob- 
tained from the difference in fitted energies of the '^Si 
(T) and ^Sq (??&) states using 3x3 fits (2x2 fits on the 
super-coarse lattices). Bootstrapping the fits together 
improved the errors in some cases, and not in others, and 
so we do not give bootstrap errors here. The results for 
the hyperfine splitting are very precise because it is the 
difference between two ground S'-wave states which are 
very well-determined. This means that statistical errors 
are small and differences between different lattice cal- 
culations show up clearly. The possibility of significant 
systematic errors in the splitting, however, must not be 
forgotten. 

Figure [T51 shows the ratio of the hyperfine splitting for 
the 1st radial excitation to that of the ground state i.e. 
{Mr' - M^'J/{Mr - M^J. The hyperfine splitting is 
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FIG. 15: Hyperfine splittings in the T and Bs systems. 
Crossed, closed and open symbols are from super-coarse, 
coarse and fine lattices respectively. Squares and triangles 
denote unquenched and quenched results respectively. The 
unquenched results are from the ensembles in Table Q] that 
include the lightest light quark mass in the quark vacuum po- 
larization at each value of the lattice spacing. Errors include 
both statistical/fitting errors from determining the splitting 
and from determining the lattice spacing scale as in Figure^] 
Additional systematic errors are 25% and discussed in the 
text. The unquenched results include a renormalisation to 
the correct value for uql (see text). The faint lines give the 
experimental limits on the — Bs splitting and the dark line 
gives the experimental result for B* — B 12^ . 



controlled by the a ■ B term in the NRQCD action and 
is then proportional to the square of C4 (see Equation^. 
C4 cancels in the ratio so in principle it can be accurately 
determined in our lattice calculation with C4 set to 1. The 
statistical errors here are large but we obtain a result of 
0.5(3). There is no sign of significant dependence on the 
lattice spacing or on the light quark mass. The result is 
in keeping with the experimental results in the charmo- 
nium system, plotted as a burst in the Figure ^Jl- The 
charmonium result is not irrelevant because we also ex- 
pect a lot of the heavy quark mass dependence to cancel 
in the ratio of splittings. 

Figure E| shows our results for the ground-state hy- 
perfine splitting in MeV as a function of the bare quark 
mass of the light quarks included in the quark vacuum 
polarisation. Wc use the lattice spacing values from the 
2S' — IS* splitting of Section IIVI to convert to physical 
units. Results are shown for the super-coarse, coarse and 
fine ensembles. The errors on the points include the very 
small statistical error from the splitting itself and also 
the much larger error from the uncertainty in the lat- 
tice spacing. This uncertainty in fact appears multiplied 
by two because the hyperfine splitting is approximately 
inversely proportional to the physical quark mass. Any 
underestimate of the inverse lattice spacing will both re- 



duce the hyperfine splitting directly and also reduce it 
indirectly by reducing the quark mass |20l |. so that the 
splitting must then be readjusted downwards to corre- 
spond to the correct meson mass. We then have a sta- 
tistical/fitting error of 5% on the ground-state hyperfine 
splitting in physical units. 

The results show significant dependence on both the 
lattice spacing (to be discussed further below) and the 
number of fiavors of quark, nf, included in the quark 
vacuum polarisation at a given lattice spacing. Neither 
of these points is surprising. In the language of a sim- 
ple potential model the hyperfine splitting arises from a 
very short-range spin-spin potential, dominated by the 
perturbative contribution which is a delta function at 
the origin. We therefore expect sensitivity to the lat- 
tice spacing and also to the value of as at short 
distance, which depends strongly on nj [s^. The re- 
sults do not however, show significant dependence on the 
bare light quark mass. This is also not surprising, as for 
spin-independent splittings, because of the large gluon 
momenta involved. Provided that the light quark mass 
is light enough, as is the case here, the splitting simply 
'counts' the number of light quark flavors. We do not 
therefore attempt any extrapolation in the light quark 
mass. 

Further examination of the dependence on the lattice 
spacing of the 2 -|- 1-flavor results is necessary, however, 
to obtain a result relevant to the real world of continu- 
ous space-time. It is important to disentangle discretisa- 
tion errors from errors resulting from unknown radiative 
corrections in C4. One important issue here is tadpole- 
improvement. We expect that radiative corrections to C4 
are not large if we use tadpole-improved gauge fields. We 
can use any sensible quantity to determine the tadpole- 
improvement factor, uq. The results will be independent 
of the quantity used in a complete calculation because 
any change will be compensated by corrections to C4. 
Here, however, we have not included any corrections to 
C4. Wc have also made slight errors in our determina- 
tion of uq from the average link in Landau gauge on the 
fine and super-coarse lattices. Because of the large num- 
ber of gauge fields appearing in the lattice discretisation 
of the B field in the a ■ B term, the hyperfine splitting 
is very sensitive to wq. It is proportional to Uq®, when 
the renormalisation of the quark mass is taken into ac- 
count |20l |. Subsequent more accurate determination of 
uql has shown that it has the value 0.8541 on the 2-f 1- 
flavor 0.0062/0.031 fine ensemble (i.e. 1% higher than 
in Table D) and 0.8102 on the 2-h 1-flavor 0.0082/0.082 
super-coarse ensemble (i.e. 0.5% lower than in Table Pi . 
We then adjust our hyperflne splitting results on these 
two ensembles by multiplying by the ratio of the sixth 
power of the mql value used in the simulation to the cor- 
rect value above, before studying the results as a function 
of lattice spacing. 

Figure [T51 shows the results for the hyperfine splitting 
on the unquenched ensembles with lightest u/d quark 
mass with this renormalisation in place. It also shows 
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again the quenched results for comparison. The discreti- 
sation effects seen are consistent with errors proportional 
to in the form: 

Hyperfinc(a) = Hypcrfinc(0)(l - {kaf + ...)• (19) 

Wc find A « 640 MeV. This variation with lattice spac- 
ing is sizeable, but smaller than previous results using an 
unimproved gluon action and an unimproved operator for 
the B field [23 • Using a constrained fitting method with 
multiple higher order terms |22l | we obtain a value for 
the hyperfine splitting in the continuum including 2-1-1 
flavors of light quarks of 61(4) MeV where the 7% er- 
ror now includes both the statistical/fitting errors on the 
individual points and the uncertainty in the continuum 
value from discretisation errors. It is dominated by the 
latter error. 

We now estimate the size of the error from un- 
known radiative corrections to C4, appropriate to tadpole- 
improvement with woL- In principle this should be done 
before studying the discretisation errors above but, as we 
shall see, the C4 effects are not dependent on the lattice 
spacing within the errors that we have so we have decou- 
pled the two effects for clarity. Radiative corrections to 
C4 are in principle 0(as) and this would give an error of 
approximately 20-30% here. However, we can make use 
of non-perturbativc information from the experimental 
spectrum of heavy-light mesons to limit this uncertainty. 

Figure 1151 also shows our results for the hyperfine 
splitting between the vector (i?*) and pseudoscalar [B^ 
ground states of the 6s system, compared to exper- 
iment The experimental results for this split- 
ting have significant errors, but results from the c-light 
mesons show very little difference between the hyperfine 
splitting for Ds and D. We therefore mark also on this 
figure the result for the B* - B sphtting (45.8(4) MeV) 
in the expectation that the B* — Bg splitting will in fact 
be very close to this. Our lattice results are obtained by 
studying correlators of NRQCD b quarks and improved 
staggered s quarks. The b quarks have exactly the same 
NRQCD action and bare quark mass as that described 
here for the T system. We used the correct value for 
UoL on the fine unquenched lattices. The bare s quark 
used for the valence s quarks (in lattice units and us- 
ing the MILC convention) is 0.04 on the coarse lattices 
and 0.031 on the fine lattices. The value on the coarse 
lattices is lower than that included in the quark vacuum 
polarization (see Table QJ, but more closely matches the 
correct value of the s quark mass determined after the 
simulation 0. 

The lattice result for the _B^*^^ — splitting is also 
sensitive to C4, but only linearly since the system has only 
one b quark. We can then use the results of Figure [T^ to 
determine C4 non-perturbatively. Using just the lattice 
Bs results we would conclude that C4=l on both coarse 
and fine lattices up to 5% errors from both the lattice 
calculation and the experimental result. We can also 
study the B* — B splitting on the lattice, but here our 
results are somewhat less precise The central value 



0.0082/0.082 



rPi - VPo 



iap2_iap; 
l^P- l^Pl 



0.018(14) 
0.016(15) 
0.005(16) 
1.1(1.3) 
0.02(1) 



TABLE XIX: P-wave fine structure for the super-coarse 2-1-1 
flavor ensemble. Energies are given in lattice units. Errors 
arc statistical/fitting only and obtained by directly from the 
fits. 



on some ensembles is as much as 10% below experiment, 
which would indicate the need for c^—l.!. On this basis 
we take a non-perturbative result for C4 of 1.0 but allow 
for a 10% error. This gives then an additional 20% error 
from C4 on the T hyperfine splitting. 

Our final result for the ground state hyperfine split- 
ting is 61(4)(12)(6) MeV where the errors are: statisti- 
cal/fitting and discretisation errors; radiative corrections 
and relativistic corrections (which we take to be 10%) re- 
spectively. Using the ratio, this gives a prediction for the 
radially excited splitting of 30(19) MeV where the domi- 
nant error is statistical. The central value of the ground- 
state splitting is larger than previous calculations using 
?!/ = 2 flavors of sea quarks |l^ I25I |2^ l27j which range 
from 20 to 50 MeV. As well as not having the correct 
value of 71/ , these older results do not have access to the 
range of lattice spacing values that wc have here or to 
the light quark masses in the quark vacuum polarisation. 
Some of them |l9l |2^ use unimproved gluon actions and 
an unimproved B fleld in the NRQCD action. These re- 
sults then have even larger discretisation errors than here 
that tend to suppress the hyp erfine splitting further [23 • 
Some of the calculations [2^ [2g do include improved B 
fields as here, and higher order spin-dependent terms in 
the NRQCD action as well. Their results indicate that 
the inclusion of the higher-order terms tend to reduce the 
hyperfine splitting [43. However, this is without any de- 
termination of C4 which is at least as large an effect. The 
new result that we give here is then more realistic than 
previous calculations and with a more reliable error. Fu- 
ture calculations need to include radiative corrections to 
C4 and relativistic corrections to spin-dependent terms. A 
more extensive analysis of discretisation effects will then 
be possible. 



B. P-wave fine structure 

Both ground and radially excited '^Po.1,2 states [xbj 
and x'bj ^^'^ ■^sll measured in the T system. No ^Pi 

(/i^,'') states have been seen. We are able to predict that 
the ground state /if, should lie within 6 MeV of the spin- 
average of the Xb states i.e. at 9900(3) (6) MeV, where 3 
is the experimental error |23| . Wc describe our analysis 
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0.01/0.05 



rPi - rPo 0.0168(20) 

l3p2 _ i3pj 0.0139(21) 



l3p_ llpj 



1.20(27) 
0.01(1) 



TABLE XX: P-wave fine structure for tfie 2+1 flavor 
0.01/0.05 coarse ensemble. Energies are given in lattice units. 
Note that these results use the half the number of configura- 
tions of other results because of a loss of data. Errors are 
statistical/fitting only and obtained by bootstrap. 







0.0062/0.031 


Quenched 


I^PI - 


-I'Po 


0.0142(34) 


0.0096(20) 


1^P2 - 


- l^'Pi 


0.0046(42) 


0.0076(20) 




^l^Po 


3.1(2.8) 


1.27(57) 


l3p_ 


-I^PI 


0.0006(20) 


-0.0002(14) 



TABLE XXI: P-wave fine structure for two of the fine lattice 
ensembles, quenched and unquenched. Energies are given in 
lattice units. Errors are statistical/fitting only and obtained 
by bootstrap. 



in more detail below, including a discussion of how well 
our results match experiment. 

The P-wave fine structure was only calculated for a 
limited selection of ensembles. For the majority of runs, 
we were unable to resolve excited splittings well and con- 
sequently show only the ground state results. The errors 
given are those from bootstrapping the fits of the individ- 
ual correlators on the coarse and fine ensembles. For the 
super-coarse ensembles they come directly from the fits. 
Tables IXIXI IXXI IXXII give the energy splittings in lattice 
units. The lattice representation for the states used 
was the E representation on the coarse and fine lattices. 
Any difference between the energy in the E representa- 
tion and in the T representation is a lattice artefact. We 
give both for the super-coarse ensemble where such arte- 
facts will be largest, but even there we cannot resolve any 
difference. The Tables also give the dimensionless ratio: 



Af (l3Pi) - M(l3Po) 
il/(l3P2) - M(l3Pi)' 



(20) 



The results obtained for this ratio have sizeable errors but 
are consistent on both unquenched and quenched ensem- 
bles with the experimental value of 1.65(9) p^ . 

Figure [TCI plots the lattice results for the ground state 
P-wave splittings, each relative to the spin average. We 
use the 25*— IS* lattice spacing values to convert to physi- 
cal units. The plot shows the coarse and fine unquenched 
and fine quenched results only since the results on the 
super-coarse ensemble have error bars which are too large 
to be useful, although they are consistent with the others. 
The Figure shows that the size of the quenched P-wave 
splittings are smaller than experiment. The results on 
the coarse 2-1-1 flavor configurations look consistent with 



> 




FIG. 16: P-wave spin splittings. Closed and open symbols 
are from coarse and fine lattices respectively. Squares and tri- 
angles denote unquenched and quenched results respectively. 
The unquenched results are from the ensembles in Table|I]that 
include the lightest light quark mass in the quark vacuum po- 
larization at the coarse and fine lattice spacing values. Errors 
are statistical/fitting only but include errors from determin- 
ing the lattice spacing as in Figure ITU There is an additional 
25% systematic error as discussed in the text. Lines give the 
experimental results where available |2^ . 



experiment, as do the fine 2-1-1 flavor results, although 
their errors are rather large. The errors shown on the 
plot are from statistics/fitting only and from determin- 
ing the lattice spacing. We double the lattice spacing 
error as for the hyperfine splitting but in general this 
is negligible compared to the statistical/fitting error of 
the P-wave splittings themselves. We now discuss the 
sources of systematic error. 

As for the hyperfine splitting, the largest sources of sys- 
tematic error potentially are the unknown radiative cor- 
rections to relevant spin-dependent terms in the NRQCD 
action. Two terms are relevant for P-wave fine struc- 
ture which complicates the picture - the a ■ B term and 
the a ■ {D X E) term - with coefficients C4 and C3 re- 
spectively. There are combinations of P-wave splittings, 
however, that can be compared to experiment to give 
non-perturbative information to fix these coefficients. 

The (J ■ {D X E) term gives rise to a coupling between 
spin and orbital angular momentum and a term in the 
mass of each state which depends on the expectation 
value of L • multiplied by C3. The a ■ B term gives 
rise to a mass term proportional to c| and either S ■ S, 
which is the same for all states, or the combination 
3(5 • nS ■ n) — S'^ which does distinguish between ^P 
states. From the expectation values of these operators in 
the different Xfe states we see that the combination 

5M(x62)-3Af(x6i)-2M(xho) (21) 
is proportional to C3 and independent of C4. The exper- 
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0.01/0.05 0.0082/0.082 

l^Dg - 1^-02 0.010(15) 0^/2*1(0) 0.5184(10) 

l^Dz - l^Dj 0.002(18) ffl3/2^^(o) 0.499(20) 



TABLE XXII: D-wave fine structure for the 2+1 flavor 
0.01/0.05 coarse ensemble. Energies are given in lattice units. 
Errors are statistical/fitting only and obtained directly from 
the fits. 

imental result for this combination is 165(4) McV and 
our lattice result on the coarse unquenched ensemble is 
164(16) MeV. This indicates that C3 is 1.0 to within an 
error of 10%. The combination 

M{xb2) - 3M(x6i) + 2M(x6o) (22) 

is proportional to c| and independent of C3. The experi- 
mental result for this combination is -46(3) MeV and our 
lattice result is -31(7) MeV. The lattice result is 30% low 
but with a large error, so that the difference is 2a. The 
limits on c| from the Bg hyperfine splitting, i.e. 1.0 with 
an error of 20%, are to be preferred for accuracy. 

We conclude from this that the systematic error from 
radiative corrections to the leading spin-dependent terms 
is 0(20%) in the ■^P fine structure. We also expect a 10% 
error from higher order relativistic corrections. Discreti- 
sation errors are not expected to be as severe here as for 
the hyperfine splitting because the P-wave states have 
a larger radial extent than the S-wave and, in poten- 
tial model language, the potentials that give rise to the 
relevant splittings are not so concentrated at the origin. 
Allowing for as much as a 10% error from discretisation 
effects, however, the total expected systematic error in 
the fine structure in this calculation is 25%. Future 
calculations will include radiative and relativistic correc- 
tions and will then resolve the '^P fine structure much 
more accurately. 

Figure ^] also shows the splitting between the spin- 
average of the ^P states and the ^Pi. This splitting is 
a result of the same S ■ S coupling that produces the 
hyperfine splitting but it is expected to vanish for P- 
wave states in a potential model because the interaction 
is focussed at the origin. We have no signal for such a 
splitting with an error of 5 MeV on the fine unquenched 
ensemble from statistics/fitting and the determination of 
the lattice spacing. Systematic errors arise from radiative 
corrections to c|, which we take to be 20% as above, and 
relativistic and discretisation errors which we take to be 
10% each. This gives a splitting of zero with an error of 
6 MeV. 



C. 73- wave fine structure 

The fine structure in the D-wave sector is on a smaller 
scale than for 5'-wave and P-wave, as expected. We 
are unable to resolve any splittings there with errors of 
around 10 MeV. D-wave fine structure was only studied 



TABLE XXIII: Results for the leading order piece of the lat- 
tice heavy-heavy vector current matrix element between the 
T and the vacuum, also known as the 'wavefunction at the 
origin'. Results are given in lattice units for the T (<I'i(0)) 
and for the T' (*I'2(0)) for the super-coarse ensembles. Errors 
are statistical/fitting only. 

on the 2+1 flavor 0.01/0.05 coarse ensemble. We looked 
at the state (using the T representation), the •^£'2 
(using the E representation) and the ^D^ (using the A 
representation) [ij. Table lXXllI gives splittings in lattice 
units with errors obtained directly from the fits. 

VIII. THE T LEPTONIC WIDTH 

The width of T decay to two leptons, Fee, can be well 
measured experimentally p3 . |45| . The matrix element 
that gives rise to the decay can also be calculated on the 
lattice and we describe that calculation here. Our most 
accurate result is for the ratio between the T' and T of 
TeeM^^^gy We obtain 0.48(5) compared to the recent 
preliminary experimental result of 0.517(10) Below 
we describe the calculation in more detail and the sources 
of error, along with improvements that will shortly be 
possible in the result. 

On the lattice Fee is calculated from the matrix el- 
ement of the bb vector current between the T and the 
vacuum. The decay width is then given by |46j : 

Fee(n5) = 16™Le?%tf^(^™a*e.)^ (23) 

where aem is the electromagnetic coupling constant, Cb 
is the charge on a 6 quark in units of the charge on an 
electron and Mx(,is) is the mass of the nth radial exci- 
tation of the T. This formula assumes a non- relativistic 
normalisation of states, (T|T) = 1. Z^atch is the renor- 
malisation constant, 1 + z^^^ofg + . . ., required to match 
the lattice vector current to a continuum renormalisa- 
tion scheme. The lattice vector current constructed out 
of NRQCD fields has leading term where tp anni- 

hilates a b quark and an anti-6. There are sub-leading 
terms which take account of relativistic and discretisation 
corrections required to match the continuum current to 
higher order. This matching can be done perturbativcly 
and is in progress The calculation is entirely anal- 
ogous to that for the B meson decay constant, whose 
Zmatch has been calculated through 1-loop 0i includ- 
ing the effect of higher order current corrections. 

For the leading order current only it is easy to ex- 
tract the matrix element from our existing results with 
no further work. It is simply the amplitude in the ??th 
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0.01/0.05 0.02/0.05 0.03/0.05 0.05/0.05 0.02 Quenched 

a^/2*i(0) 0.30149(66) 0.29830(59) 0.29441(67) 0.29142(78) 0.28563(53) 0.26152(50) 
a3/2^a(0) 0.2446(66) 0.2401(50) 0.224(14) 0.237(11) 0.230(11) 0.2346(48) 



TABLE XXIV: Results for the leading order piece of the lattice heavy-heavy vector current matrix element between the T and 
the vacuum, also known as the 'wavefunction at the origin'. Results are given in lattice units for the T (5'i(0)) and for the T' 
(^'2(0)) for the coarse ensembles. Errors are statistical/fitting only. 



0.0062/0.031 0.0124/0.031 Quenched 



a3/2*i(0) 0.17666(49) 
0^/2^3(0) 0.1357(24) 



0.17796(44) 0.15159(31) 
0.1372(19) 0.1298(17) 



TABLE XXV: Results for the leading order piece of the lattice 
heavy-heavy vector current matrix element between the T and 
the vacuum, also known as the 'wavefunction at the origin'. 
Results are given in lattice units for the T (*I'i(0)) and for the 
T' (^'2(0)) for the fine ensembles. Errors are statistical/fitting 
only. 




0.02 0.04 0.06 0.08 0.1 
miig,,(GeV) 

FIG. 17: The ratio of leptonic width times mass squared for 
the 1st excited state of the T to that of the ground state. 
Crossed squares are from super-coarse unquenched ensem- 
bles, closed squares from coarse unquenched ensembles and 
open squares from the fine unquenched ensembles. The closed 
and open triangles are from the coarse and fine quenched en- 
sembles respectively. The closed diamond is from the coarse 
Uf = 2 ensemble. The burst represents the current experi- 
mental result. 



T state of the local (delta) smearing, i.e. a{loc,n) from 
equationlHl appropriately normalised. In potential model 
language this is known as the wavefunction at the origin, 
^'(0), and is obtained from the lattice results when the 
correlators are divided by the product of the number of 
colors (3) and the number of non-relativistic spins (2). 
Then \l/(0) — (T„| J„|0)/-\/6, accounting for the factor of 
6 in the denominator of Equation [53 Results for this 
quantity are given in lattice units in tables IXXIIII IXXIVI 



and IXXVI for the super-coarse, coarse and fine ensembles 
respectively. 

There will be an 0{as) systematic error from Z match 
if we try to calculate Te+e- directly from these numbers. 
Instead we form the ratio: 



ree(25)Af.2 



T(2S) 



ree(l^)M2 



T(1S) 



1^2(0)1^ 

l*i(0)|^ 



(24) 



in which the Zmatch factor for the current cancels. We 
then expect to be able to determine this ratio from the 
leading current with errors coming from the matrix ele- 
ments of relativistic corrections to the current. 

Figure El shows this ratio for our lattice results on 
super-coarse, coarse and fine ensembles. The current ex- 
perimental result |2^ is given by the burst. The lattice 
results in the quenched approximation arc in strong dis- 
agreement with experiment. Previous quenched results 
also indicated this, but less precisely |l3.l20| . The results 
move down towards experiment once light quark vacuum 
polarization effects are included and the light quark mass 
falls. However, the results on the unquenched ensembles 
with lightest light quark mass are still some way from ex- 
periment. The plot shows no dependence on light quark 
mass for the most chiral ensembles, as expected. We 
therefore do not attempt a chiral extrapolation but now 
look at lattice spacing dependence for those results with 
lightest light quark mass. 

The quenched results for the ratio and those from the 
most chiral unquenched ensembles are plotted against the 
square of the lattice spacing in Figure^l This plot shows 
a significant dependence on a^, indicating discretisation 
errors in the results. These may come from the current 
operator or from the action at the level of asO^ and it 
is not clear which is the most important at this stage. 
The wavefunction at the origin is a very short-distance 
quantity (akin to the hyperfine splitting) and it is reason- 
able for it to show significant discretisation effects when 
other quantities, such as radial and orbital splittings in 
the spectrum, do not. 

We can directly calculate the leptonic width for the T 
and T' in the absence of Zmatch to see how each behaves. 
We use the formula given in Eauation l23l taking OLemi'mb) 
= 1/132 m. We obtain a value ree=1.43(7) keV for the 
T from the super-coarse 2-|-l flavor unquenched ensem- 
ble, 1.32(1) keV from the coarse 2+\ flavor 0.01/0.05 
unquenched configurations and 1.28(1) keV from the fine 
2+1 fiavor 0.0062/0.031 configurations. The error given 
here comes only from the error in the determination of 
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the (25' — IS) lattice spacing which appears cubed in 
this result. The fine lattice numbers is slightly lower 
than the recent preliminary experimental result, 1.34(2) 
keV The absence of Zmatch or higher order correc- 
tions gives a sizeable additional systematic error to the 
lattice result, but there is at least nothing to indicate 
that large systematic shifts arc required to match exper- 
iment. The quenched results for the T arc significantly 
lower at 1.0 keV on the fine lattices. For the T' there 
is more sign of discretisation in the bare leading-order 
matrix element. From this, we obtain Fe+e- = 1-2(1) 
keV on the super-coarse unquenched ensemble, 0.77(5) 
keV on the coarse unquenched lattices and 0.67(3) keV 
on the fine unquenched lattices, with errors as above. 
This is not inconsistent with the fact that the excited 
state has more structure in its 'wavcfunction' and there- 
fore might be more susceptible to short distance errors 
than the ground state. The fine lattice result is slightly 
higher than the recent preliminary experimental value of 
0.616(13) keV HI. 




(GeV") 



FIG. 18: The dependence on a? of tlie ratio of leptonic width 
times mass squared for the 1st excited state of the T to that 
of the ground state. Crossed, closed and open squares are 
from the most chiral of the super-coarse, coarse and fine un- 
quenched lattices respectively and closed and open triangles 
are from quenched lattices. The burst represents the current 
experimental result, and the cross below this is our contin- 
uum extrapolated result, including an estimate of the leading 
relativistic corrections as described in the text. 

As for the hyperfine splitting earlier (Equation we 
consider a-dependence in the ratio of Equation |31 of the 
form 



ratio(a) = ratio(0)(l + {Kaf 



(25) 



However, here the a-dependence is stronger than for the 
hyperfine splitting (A « 900 MeV) and the super-coarse 
results are therefore correspondingly less reliable. Using 
again a constrained fit we find a continuum value for the 



ratio of 0.50(4) for the leading order current. The 8% 
error is a combination of the statistical/fitting error and 
the effect of discretisation errors. 

The higher order relativistic current corrections are 
given by operators of the form [50l | : 



1 



8(M0)2 



(A(2)tcr -I- (7A(2) - 2(cr • At)fT((T • A))V'. (26) 



For an T at zero momentum the A'^' operator acting 
on cither quark or anti-quark will give the same matrix 
clement. This is similarly true for the operator made 
from one derivative on each. The zth component of the 
correction operator then becomes: 



1 



4(M0)2 



(-A(2)t(j^ + 2Aj'T,)V'. 



(27) 



where Dij = AiAj — A^^^/S. The second term above is 
the piece which couples to a D-wave state Q and gives 
a leptonic width to the ^Di state. The first term is the 
one that we are concerned with and that gives a relativis- 
tic correction to the S'-wave leptonic width of the form 
x1'A(2)cr/6(M0)27/;. We calculate the matrix clement of 
this term by inserting a A*^^^ operator at the correlator 
sink. We then fit the correlator with the insertion si- 
multaneously with the other correlators used to get the 
matrix element of the leading order current above. We 
have done this for the 0.01/0.05 2-f 1 fiavor coarse ensem- 
ble only. There the current correction matrix element (in 
lattice units) is -0.0107(2) for the T and -0.0131(4) for the 
T', i.e. a 3.5% and 5.4% correction to the leading matrix 
element respectively. This induces a further 4% reduc- 
tion in the ratio of leptonic width times mass squared 
for the 2S to IS*. The size of the relativistic corrections 
and the overall shift from their difference is perfectly in 
keeping with expectations based on power-counting in . 
Higher order relativistic corrections are then likely to be 
negligible. 

We apply this shift of 4% to our continuum leading 
order current result. Wc also add an additional 4% error 
(a^w^) to allow for missing radiative corrections to the 
relativistic corrections in the ratio (radiative corrections 
to the leading order term cancel, as discussed earlier). 
This gives our final answer for the ratio of Equation [21 
of 0.48(5), to be compared to the recent preliminary ex- 
perimental result of 0.517(10). Both of these are also 
marked on Figure [T51 Our result is compatible with ex- 
periment but not as precise. Our 10% error is dominated 
by the fact that the lattice results depend strongly on the 
lattice spacing. 

The lattice calculation of Fee can be significantly im- 
proved. Once the renormalisation of the leading and sub- 
leading lattice NRQCD current operators is calculated, 
we will be able to compute the Fee results for each ra- 
dial excitation of the T separately with systematic er- 
rors at the 10% level coming from unknown terms in 
Zmatch E3- We will also discover whether the current 
corrections contain a sizeable contribution in the form of 
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FIG. 19: The ratio of lattice results to experiment for var- 
ious mass splittings involving b quarks. The results are 
for quenched simulations (left) and simulations incorporat- 
ing 2-1-1 flavors of sea quarks (right). All the unquenched 
results come from this work. We have taken results from 
the most chiral coarse lattice ensemble since those have the 
best errors in general. Errors include statistical/fitting er- 
rors for the quantity in question and those from determin- 
ing the lattice spacing. The top three quenched results come 
from this work, but the lower two are estimates from earlier 
work. The quenched ID — IS splitting is taken from |53 | 
using the state. The quenched 2ms^ — mr splitting is 
estimated from fs^ . adjusted to use the lattice spacing from 
the T 25 — IS spacing. 

a discretisation correction, so that the scaling with lattice 
spacing improves. As discussed earlier, the ratio of lep- 
tonic width times mass squared will be more precise than 
this because the errors will be set by corrections to the 
current and errors in the current correction piece. 
These errors should be at the level of a few percent. 

IX. CONCLUSION 

We have presented a state-of-the-art calculation for T 
spectroscopy in lattice QCD. The major improvement 
over previous calculations is that we now have results on 
gluon field configurations that include 2-1-1 flavors of sea 
quarks with the u/d quark masses within reach of their 
physical values. This means, as we demonstrate here, 
that it is now possible to fix the parameters of the QCD 
Lagrangian unambiguously. The T system is a good 
place to test this because of the large number of well- 
determined excited states. Figure IH^ shows a 'ratio' plot 
of lattice results from this paper compared to experiment. 
The unquenched results on the right are in striking agree- 
ment with experiment compared to the quenched results 
on the left. This ratio plot differs slightly from that in Q 
for the quenched results because of improved fits for the 
quenched ensembles here. 



Our results enable us to determine physical values for 
the parameters ro (0.469(7) fm) and n (0.321(5) fm), 
often used to compare lattice spacing values between dif- 
ferent ensembles in lattice calculations. 

Further improvements to the calculation of radial and 
orbital excitation energies will require further improve- 
ments to the NRQCD action, in particular the calcula- 
tion of radiative corrections to v'^ and terms. This is 
in progress. Statistical precision can also be improved 
since double the number of configurations used here are 
available as well as further configurations with lighter sea 
quark masses 0, . Precision on P and D-wavc states 
would be improved by better smearing. 

Our result for the b quark mass in the MS scheme 
at its own scale is 4.4(3) GeV, using a 1-loop matching 
of the lattice to the continuum. A 2-loop calculation of 
the heavy quark self-energy would reduce the error signif- 
icantly on the b quark mass as well as make the remaining 
error estimate more reliable by fixing the scale for in 
the matching coefficient. This is currently underway. 

The fine structure in the spectrum is only calculated to 
leading order here, so sizeable systematic errors remain. 
We arc able nevertheless to predict the ground-state hy- 
perfine splitting more accurately than in the past and the 
radially excited hyperfine splitting for the first time. Our 
results are 61(14) MeV and 30(19) MeV respectively. To 
reduce the errors requires, as above, the determination of 
radiative corrections to the leading spin-dependent terms 
in the NRQCD action and the inclusion of sub-leading 
terms. Further improvement to reduce discretisation er- 
rors is also an issue. 

The calculation of the Icptonic width for T and T' is an 
important example of a simple decay rate that can be cal- 
culated in lattice QCD, without the ambiguity that arises 
from potential models [sJl- Our results on unquenched 
configurations for the ratio of leptonic widths (multiplied 
by squared masses) for T' and T show a much closer 
agreement with experiment than those on quenched con- 
figurations. Sizeable discretisation errors prevent us giv- 
ing a very accurate result here, but wc obtain 0.48(5) for 
the ratio. Further precision will come from the calcula- 
tion of the renormalisation constant required to match 
the lattice and continuum currents responsible for the 
decay. The calculation will yield an improved current 
operator with reduced discretisation errors. This is also 
underway. 
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